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Abstract. The magnetic extension of the Thomas-Fermi-Weizsacker kinetic energy is used within density- 
functional-theory to numerically obtain the ground state densities and energies of two-dimensional quantum 
dots. The results are thoroughly compared with the microscopic Kohn-Sham ones in order to assess the 
validity of the semiclassical method. Circular as well as deformed systems are considered. 
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1 Introduction model similar to that developed by Zaremba and cowork- 



ers |^| . It is our aim in this paper to extend those calcula- 



I 



^ . Semiclassical approaches to many-body systems are a very tions by including the effect of a perpendicular magnetic 

c : 

O valuable tool since they provide physical insights which ne id B using the magnetic extension of the kinetic energy 

j>! ; otherwise are very difficult to achieve. In fact, they have within density-functional theory. As in Ref. H] we will pay 
•i—i . 

X ; been applied since many years ago to describe different specia i attent ion to the quantitative comparison with the 

— ■ systems such as atoms, atomic nuclei, metals and, more microscopic Kohn-Sham solution in order to asses the ac- 

recently, metallic clusters and electronic nanostructures. curacy limits of the TFW densities and energies as a func- 

Two dimensional quantum dots are not an exception and t ion of B. The magnetic extension of the Thomas- Fermi 

have been analyzed using the Thomas-Fermi models in, theory was rigorously presented by Lieb et al. § and a 

for instance, Refs. (I],|^,|3|,[|,[|,|(|. selfconsistent numerical application to circular dots was 

In Ref. we performed calculations for quantum dots done in Ref. but, to our knowledge, no detailed corn- 
using a selfconsistent Thomas-Fermi-Weizsacker (TFW) 
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parison with microscopic calculations has been given in 2 Magnetic TFW functional 

the literature. We also present in this manuscript symme- 

Using a local approximation to density functional theory 

try unrestricted calculations for deformed dots that had 

we assume that the energy of the system can be written in 

not been considered before within this model. 

terms of the electronic spin densities /Ojj(r), where n =f, j, 

A peculiar characteristic of the magnetic Thomas-Fermi as e = J dr£[p hPl ;B}. Notice also the explicit depen- 

functional is its first-derivative discontinuity in the density dence on the magnetic field B which is considered as a 

dependence ||, at particular density values. The physical functional parameter. The different contributions to the 

origin for this is found in the formation of constant en- energy density may be written as 

ergy Landau-bands in the non-interacting Fermi gas at ^ 

£ [p hPl ;B}^T[p r ,B}+T[p l ;B} + -v H (r)p 

certain magnetic fields. This implies that the mean field is 

,. - ( .v t t . + £xc(PhPl) + v C xt(r)p + £z(PhPi; B ) , (!) 
discontinuous, thus manitcstmg trie formation ot incom- 
pressible regions in the system, each one characterized by where the nrst two pieces give the k ine tic energies of spin 
the number of full Landau bands, i.e., the filling factor up and down electrons; the third is the Hartree energy in 
v. It is worth to point out that experimental evidences terms f the Hartree potential v H and the total density 
of incompressible stripes at the edges of quantum dots p = P] + Pi; the fourth is the exchange-correlation con- 
and antidots have been obtained by means of far-infrared tribution; the fifth is the energy due to the external po- 
spectroscopy g. As we will show, the semiclassical Euler- tential ?w and, the last one corresponds to the Zeeman 
Lagrange equations selfconsistently determine the density contribution. As in Refs. §,§, the kinetic energy contains 
profile, energy and chemical potential of the quantum dot, a pure Thomas-Fermi term and a gradient (Weizsacker) 
which turn out to be in an overall good agreement with one5 T = TTF+Tw _ The gradient term is given by 



full Kohn-Sham results for increasing magnetic fields up fe 2 fX7 \2 

to v = 1. 2m Pn 



MP,} = ^A^^ , (2) 



Section 2 of the paper is devoted to the presentation 



with A = 1/4, while the Thomas-Fermi piece is Bf 



of the energy functional as well as the minimization equa- ttf[Pt)', B] = — hui c D 

tions. In Sec. 3 the results for circular as well as deformed _|_ ^ j g _|_ £ j s / ) i ( ! i 

quantum dots are given and, finally, Sec. 4 presents the 
conclusions. 



In this last expression uj c — ^| is the cyclotron frequency, 
D = 2^h c lSi ^ ne Landau level degeneracy per unit area 
and S v — [^j] is the integer part [|| of the local filling 
factor v-q — and gives the index of the highest fully 
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occupied Landau band. The two contributions in Eq. (^|) equations read 

give therefore the kinetic energy of the fully occupied Lan- ( se ( se 

Sp ~ V J <5pt ^ 

dau bands and that corresponding to the last partially \ = I ' 



SE. _ n \ SE_ 

Sm — U ^ S Pl " A 4 

For convenience, we choose to work with the second set 
2 . .p-i which can be transformed, introducing new variables ip v — 



filled band, respectively. It can be shown |3|] that, in the 
limit B — ► 0, the non-magnetic Thomas-Fermi functional 
ttf{p, B = 0) = ^i^p 2 is recovered from Eq. 

For the exchange-correlation energy £xc we have used 

the LSDA functional based on the Tanatar-Ceperley cal- _ 4A 2^ V ^ + I Wext + VH+ 
culations for the 2D uniform electron gas ju| and the von + ^ + (Xi ^L. [1b b^ ^ = M ^ , (6) 

Barth-Hedin interpolation for intermediate polarizations 
pj| . The expression can be found, e.g., in Ref. [Q, where 
Kohn-Sham results for parabolic dots at = were given. 



p v , into the alternative Schrodinger-like equations 



where we have defined a-\ = 1, ai = — 1 and also intro- 
duced the contribution from the Thomas-Fermi energy 



Notice that current-density dependence is not included in C v — hui c ( S v + — ) (7) 



2 



the functional. This could in principle be done within the 
so-called current-density-functional theory although it is 
known that the contribution to the ground state energy 
from these terms is in general quite small and only at very 
high magnetic fields they can be of relevance [0,[Ilt]. 



The solution of the two coupled Eqs. (||), i.e., for each 
spin component, has been obtained numerically by dis- 
cretizing the two dimensional xy plane into a uniform grid 
of points and using the imaginary time-step method. The 
grid size is typically 70 x 70 points, while the Laplacian 
The Zeeman energy reads, in terms of the effective operator is discret i Z ed by using 7 points formulas. The 
gyromagnetic factor g* and Bohr magneton (i B , stability of the results when increasing these values has 

been ckecked. The Kohn-Sham results, we will compare 
£z{p-\ , Pi', B) = — g* [is B (pf — pi) . (4) with below, have been obtained using a similar method 

developed by us in Ref. 

The energy functional is minimized by the ground state 
spin densities, or equivalently by the ground state total ^ Results 
density p and magnetization m = p\ — pi , with the con- 

3.1 The TF plateaus 

straint of conservation of the total number of particles. 

The corresponding Lagrange parameter p is, by defini- We begin the results section by discussing the effect of the 
tion, the chemical potential. The two sets of equivalent discontinuous contribution to the potential. We present 
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here calculations for two different dots: a circular one con- down the second band is only partly occuppied. When go- 

taining N = 42 electrons under parabolic confinement ing towards the edge (increasing £) the second spin down 

^oxt = 1/2 uj^r 2 |l6|,|l7j and a second one containing N = band is rapidly depleted and at a larger r the same hap- 

20 electrons in a deformed parabola pens with the second spin up band. This behaviour of the 



/ \ 1 2 4 / 2 . a i 2\ / \ microscopic solution is in excellent agreement with that 

"oxt(r) = ^ + ma + P V ) , (8) P 6 

inferred from the plateaus of the left panel, thereby show- 

with anisotropy factor (3 = 0.75 and a coefficient ujq 

ing the quality of the model. The formation of the plateaus 

given by N p = 20 and the same r s as for N — 42. 

is not as clear in the KS densities because of rather large 

The spin up and down densities for the circular dot are 

density oscillations, quite similar to the Friedel oscillations 

displayed in the upper-left panel of Fig. [l], in comparison 

found in metals. 

with the corresponding Kohn-Sham ones. The horizontal 

The two lower plots of Fig. |l| show the plateaus in spin 

dotted lines indicate the values p = D and 2D where D 

up and down densities for the deformed dot at B = 1 T. 

is the density for bulk filling factor v = 1 (Sec. 2). The 

A behaviour similar to that of the circular case is inferred, 

TFW density is clearly giving flat regions, i.e., plateaus, 

although in this case the plateaus adjust their shape to the 

at the densities of the bulk integer filling factors, which 

anisotropy j3 = 0.75 of the confining potential. At larger 

must be associated with incompressible stripes in the fi- 

magnetic fields, however, deviations are obtained as we 

nite system. Due to the effect of the Weizsacker term the 

will show when presenting the systematic results in Sec. 

transition between different plateaus is quite smooth. In 

3.3 

this particular case the spin down density has attained the 
first plateau and is beginning to fill a second one around 

3.2 Circular dot with 42 electrons 

the dot center. On the other hand, the spin up density 

has already reached the second plateau at the center of In this subsection we show in a systematic way the re- 
the dot. suits for the dot containing N = 42 electrons in a cir- 

The correlation of the plateaus with Landau bands is cular parabola, with r s = 1.5 ag and N p = 42. Figure 
made even clearer in the upper-right panel, where the || displays the evolution with magnetic field of the den- 
Kohn-Sham eigenvalues are plotted as a function of or- sity and magnetization profiles in comparison with the 
bital angular momentum £. In this plot the horizontal line Kohn-Sham ones. In general the TFW density and mag- 
at e w — 1.24H* indicates the Fermi energy. A proportion- netization are correctly averaging the KS values with a 
ality between £ and r may be established by noting that rather good agreement for all the magnetic fields consid- 
high i values imply outer orbits. Therefore, at the dot ered. At B = 5 T the TFW correctly yields equal density 
center (low ^'s) two spin up bands are filled while for spin and magnetization distributions, due to the achievement 
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of full polarization. At this magnetic field the KS result Figures g and [f| show within TFW and KS, respec- 
corresponds to the maximum-density-droplet (MDD) so- tively, the local filling factors v for selected values of the 
lution p9| , in which the single-particle angular momenta magnetic field. A similar behaviour as that discussed for 
are succesively occupied up to the ^ max — N —I value. In- the circular case is obtained. In fact we recognize the pref- 
creasing the magnetic field still further, i.e., entering the erence of the TFW density to produce flat regions asso- 
region of fractional filling factor v < 1, the dot evolves ciated with the integer filling factors. By looking at the 
by reconstructing the edge, as seen in the B = 6 and 7 T central region we can also identify the progressive deple- 
panels. However, this physical behaviour is apparently not tion of the central plateau when increasing the magnetic 
well reproduced by the TFW model. field. For instance, one can follow the evolution v = 3, 2, 1 
In the upper-left panel of Fig. | a quantitative com- for B = 1,1.5,2.5 T, respectively. This behavior is also 
parison between the energy per particle in TFW and KS inferred from the KS results (Fig. |) although it is some- 
is provided. We see that the difference remains below 2%, how masked by the large oscillations of quantum origin, 
although on the plot it may seem magnified because of the A s in the circular case, the prediction of the dot polariza- 
expanded scale. In the upper-right panel a comparison of tion with magnetic field is also in good agreement with 
the TFW chemical potential n with the KS Fermi energy the microscopic result. This can be qualitatively seen in 
is given, which is again indicating the good estimate given Figs, g and | by noting the clear depletion of the spin 
by the TFW of the less bound electron, until the edge re- down density starting at B = 3 T which ends with a fully 
construction begins. We mention that we remain here in polarized dot (S — 10) for B > 4.5 T. 
the limit T — > 0, although a small value of T is sometimes 

A conspicuous prediction of the KS model is the grad- 

necessary to converge the KS results. 

ual change in shape of the quantum dot when increasing 
the magnetic field. This is quite evident for B > 4 T, 

3.3 Elliptic dot with A = 20 when the dot can no longer be considered elliptic, but 

rather rectangular in shape. We attribute this to the com- 

The lower panels in Fig. ^ represent the comparison of the petition between the symmetry of the external field and 

ground state energy per particle E/N and TFW chemical the preference for circular edge reconstructions induced by 

potential vs KS Fermi energy corresponding to the 20- large magnetic fields. A more abrupt transition to a cir- 

electron dot in a deformed parabola (Sec. 3.1). As for the cular shape was obtained in Ref. [2(J within the ultimate 

previous circular system the agreement between the TFW jellium model, which permits the deformation of the ex- 

and KS ground state energies is remarkable, not exceeding ternal potential. In our model this is fixed and, therefore, 

in this case a 3%. it seems natural that a stronger competition is present. 



6 LI. Serra, A. Puente: Magnetic Thomas-Fermi- Weizsacker model for quantum dots . . . 

When comparing with the TFW results, we notice that this model is not able to yield the shape changes induced 
although the central plateau indeed seems to evolve to- by the magnetic field, as found in the KS result. 



wards a rectangle (for B = 4.5 — 5.5 T), the outer edge 
remains always elliptic. The defhciency of the semiclassical 
model in reproducing morphological changes attributed to 
the magnetic field can be explained by the fact that within 
TFW the magnetic field effects are taken into account in a References 
purely local way (see Eq. (0)), and thus can hardly induce 
any influence on the global shape. Another missing feature 
in the TFW results is the incipient electron localization 
seen in the KS panels for B > 6 T (Fig. 5). This localiza- 
tion has also been predicted within Hartree-Fock theory 
|plf and more recently, using current-density-functional 
theory fTl. 



This work has been performed under Grant No. PB98- 
0124 from DGESeIC, Spain. 
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